TUM-T3 1-43/93 
The Decay B K*-f from QCD Sum Rules 



Patricia Ball 

Physik-Department, TU Miinchen, D-85747 Garching, FRG 
(August 7, 1993) 



Abstract 

Using the QCD sum rules approach and the complete leading-logarithmic 
short-distance coefficients determining the penguin-amplitude, we obtain 
B{B K*-f) = (6.8 ±2.4) • 10"^ for top-quark masses in the range (150-200) 
GeV and \Vcb\ = 0.035 for tb = 1.5 ps. The ratio to the inclusive rate is 
B{B ^ K*-f)/B{B Xs-f) = (20 ± 6)%. 



Since some time, the penguin-induced decay b sj has received continuous interest 
as a possible probe of new physics (cf. e.g. [0-0]). The interest quite naturally increased 
when experimental results on the exclusive radiative decay B — > K*'y became available 
The interpretation of these results, however, is not obvious, since existing calcula- 
tions focused mainly on the inclusive decay B — > Xs'j with inidentified strange hadrons 
Xs in the final state. The rate of this decay can be obtained from the perturbatively cal- 
culable short-distance expansion of the effective Hamiltonian. The existing estimates of 
hadronic corrections rely mainly on quark model [§-0] as well as QCD sum rule calculations 
PJ^, and yield rather different results. To be specific, theory is asked to provide one with 
the hadronic matrix element ( K*{pk*) \ sa^^q^Rh \ B{pb) ) at vanishing momentum transfer 
squared, = {pb ~Pk*Y = 0, where R = {1 +75)72 is the usual projector on right-handed 
quarks, and a^u = i/2 [7^,7i,]_. In this letter, we aim at to give an updated analysis of the 
exclusive decay B —>■ K*^, using both short-distance coefficients in leading-order accuracy 
im and improved QCD sum rules.Q We will not explain the QCD sum rules approach in 



detail, but just mention that it has established itself as a reliable tool to infer the gross fea- 
tures of QCD induced non-perturbative dynamics of hadronic matrix elements. The method 
relies on the field-theoretical aspects and features of QCD and was designed to make max- 
imum use of known manifestations of non-perturbative QCD. Originally invented for the 
calculation of vacuum-to-meson transition amplitudes |11| , it soon found application to the 



calculation of the electromagnetic form factor of the pion |12| and other meson-to-meson 



^ After finishing our calculations, we became aware of Ref. |2^ where it is claimed that the Wilson 
coefficients calculated in |1C] are not completely correct. The change in numerics is, however, 
negligible (~ 1%). 
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transition amplitudes (cf. for a review). Although QCD sum rules in general yield less 
detailed results than fine-tuned models, they have the advantage that only a small number 
of parameters is needed that have an evident physical meaning (e.g. quark masses) and/or 
characterize the non-perturbative regime of QCD (e.g. the so-called quark condensate, the 
order parameter of chiral symmetry breaking). Once these parameters are fixed from well 
known processes, they can be used to calculate for instance heavy meson decays. 
The decay we are interested in is determined by the effective Hamiltonian 

= ^v,bv:,Y.c^{f^mfi), (1) 

i=l 

describing an effective five-quark theory where the effects of the top-quark and the W-boson 
are integrated out. The operators Oi mix under renormalization, their anomalous dimen- 
sions are known completely to leading-order and partly to next-to-leading order accuracy 
||l| J10|JT^ . Since a calculation of 0{as) corrections to our sum rules is beyond the scope of 



this letter, we for consistency only use the leading-order coefficients calculated in ITO 
The only operator relevant for the decay B — * if* 7 is 

07 = sa,,{m,R + m,L)bF^''. (2) 

Here rrib and are the masses of the b- and s-quark field, respectively, F^i, is the elec- 
tromagnetic field strength tensor, L = (1 — 75)/2 the projector on left-handed quarks. In 
calculating the decay rate, we note that the structure (7,^,^75 is not independent of a^u, 
but can be expressed as cx^vTs = —i^fiuafBO'"^ /2, and thus the form factor decomposition of 
( K*, A I sa^^q^h \ B ) determines the decay rate completely. With 

( if *, A I m^.q^'h \B) = -iF{q\ fi) e^.^^eSJ^^^., (3) 

where q = Pb — Pk* is the photon momentum, A denotes the helicity state of the K* and 
e^^^ is the corresponding polarization vector, we find 

r(5 - if*7) = T7^,Gl \V:M'C'rif^)F\0,^,){ml + ml) ^"^^ ~ . (4) 

In (^ and @, we have made explicit the scale-dependence of the matrix element by ascribing 
the form factor a dependence on the normalization scale fi. The inclusive rate is given by 

cy 

TiB ^ X,7) = i^,Gl |V;;K,P C'.ifi) ml ml, (5) 

neglecting the s-quark mass. 

In the framework of QCD sum rules, the form factor F{0, fi) can be calculated as ground 
state contribution to the three-point correlation function 

Tau{p%,PK*^(l^ = 0) = 

= t^[d'x d%j e-'^^^^'P'<'y ( I T q{y)ia.s{y) s(0)g'^a^,6(0) h{x)i^^q{x) \ ) 

= iT{pl,pl,)eau(3^p^BPh- (6) 
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T can be expressed as (double) dispersion relation where the double spectral function is 
given by the contribution of intermediate on-shell particles coupling to the currents. In 
zero-width approximation we find 

Tau = -I F(0, /i) 1^ -2 -2^T-2 + contmuum, (7) 

A \Pb " ^b)\Pk* ~ ^K*) 

where the first term on the right-hand side just contains the expression we are interested in, 
the second one stands for the contributions of higher resonances and many-particle states. 
Introducing the meson leptonic decay constants Jb and fx* for the B and the K* meson, 
respectively, as 

{ I gz75& \B) = ( I q^,s \ K\\) = mK*fK*e\ (8) 



rrib 



we find 



rjn ^ T^K*fK*m%fB 1 , . 

1 = -T (U, /i) —T, TTT^ ^r~: + contmuum. (91 

rUb [p% - m% ) [pj^, - mj., ) 

Using quark-hadron duality, we model the continuum part by the purely perturbative double 
spectral function above certain thresholds in the dispersion variables, the so-called contin- 
uum thresholds, s% and , respectively. In order to enhance the contribution of the ground 
state and to suppress the continuum contribution, we subject T to a Borel transformation 
B (cf. Inl) yielding 

(1°) 

By means of this transformation, the variable gets effectively replaced by M^, the Borel 
parameter. Evidently, large values of the dispersion variable s, i.e. the continuum contribu- 
tion, get exponentially suppressed. 

Whereas T is calculable within pure perturbation theory for p^, — 00 only, the 

QCD sum rules approach allows one to get closer to the more interesting physical region 

— ml ^ p'j^t — ml ~ 0. Performing an operator product expansion of (^, it proves 
feasible to parametrize the unknown long-distance behaviour by vacuum expectation values 
of certain gauge-invariant operators, the so-called condensates . These matrix elements 



vanish, when taken over the perturbative vacuum, but acquire finite values over the physical 
vacuum. Thus T can be expressed as 

npWK^)=Y.T^''\vWK*){On). (11) 

n 

The coefficient functions T*^") contain the short-distance behaviour of the correlation func- 
tion and are calculable perturbatively. In this letter we take into account all operators up to 
dimension n=6, i.e. the quark, gluon, mixed, and four-quark condensates. The evaluation of 
these contributions proceeds along standard lines, for the coefficient of the gluon condensate 



we use the method proposed in |15]. All coefficients are calculated to O(a^) and including 



their full dependence on m^. For lack of space, we do not give the formula in this letter. 
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Applying the Borel transformation in both the variables and pj^, , (j^) yields the final 
sum rule 



F(0,/i) = J^' gm|/M|+<./M|, MlMl.f{Ml,Ml,). (12) 

mK*jK*m%jB 

Here T is the correlation function with subtracted continuum contribution, and Mk* are 
the Borel parameters. At this point, a few remarks are in order. First, the quark mass rrib 
appearing in Eq. (0) is the renormalization scheme and group invariant pole mass, related 
to the running mass in the MS scheme by 

m„„,^™W.)(l + ^{5 + >n4:)). (13) 

In the numerical evaluation we use the value 771^ = 4.8 GeV. For the mass of the s-quark, 
we use the running MS mass with rris^fi = 1 GeV) = (0.20 ± 0.05) GeV. Actually, the sum 
rule proves not very sensitive to that value. For the condensates, we use the standard set of 
values (at the renormalization point /i = 1 GeV): 

{qq) = (-0.24GeV)3, {qga^^^^'^q) = 0.8 GeV^qq), 

16 

{aJ'K'G^) = 0.012 GeV\ 7r«.(g7"A^g Yl ^7rA^g) ^ -y v^«.(gg)^ 

u,d,s 

47ias{duud) ^ Anasiqqy. (14) 



For fx*, we use the experimental value extracted from r K* u^, fx* = 0.21 GeV [16 



For Jb, we use a two-point sum rule (see [|r^, e.g.) without radiative corrections. In the 
limit of infinitely heavy quarks, this procedure is required to obtain the right normalization 



of the form factor [|T8l, and largely reduces the dependence of the sum rule on the Borel 



parameter and the infiuence of (in our case unknown) radiative corrections [|T^]. We take over 
the procedure to the present case whith one light quark, and use half the value of the Borel 
parameter in the sum rule for Jb and equal continuum thresholds s% in both the two- 
and the three-point sum rule. In addition, we evaluate (|12]) at a fixed ratio M'^/M]^* = 3. 
The numerical values of the continuum thresholds are taken from the region of maximum 
stability in the two-point sum rules, where one finds s^^. = 1.7 GeV^ and s% = 34GeV^, 
and a "sum rule window" TGeV^ < M| < 10 GeV^ (cf. |2gj2T|). 

Next, we have to say some words about the residual scale-dependence of F(0,/i). In 
contrast to calculations at partonic level, the condensates induce a scale-dependence of the 
sum rule already at leading order in as- Whereas previous sum rule calculations tacitly 
assumed fi = rrih we argue that the effective scale is much below that. This becomes 
plausible when one takes into account that the sum rule calculation is an intrinsic off- 
shell one where the correlation function is calculated in the not so deep Euclidean region 
of external momenta and then analytically continued to the physical region. Thus the 
characteristic scale is rather given by the virtualities of the particle currents than by a 
singularity of the correlation function at threshold. Our reasoning gets support from the 
analysis of the heavy quark limit where the use of a low normalization scale M^/(2mQ) ^ 
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1 GeV in the radiative corrections to the sum rule for fs is absolutely stringent [^]. From 



all that we conclude that the proper normalization scale of the sum rule is coupled to the 
Borel parameters. In view of the successes of previous form factor determinations ( ||20|| , 
e.g.), we conform to the choice made there and use n = {M'^Mj^*)^/^ which typically is about 
2 GeV. Note that a further discussion of that point would require the complete calculation of 
0{as) corrections to the correlation function which is complicated by the effects of operator 
mixing. 

The next step now is to disentangle the intrinsic M^-dependence of the sum rule for 
the form factor reflecting the inherent uncertainty of the QCD sum rule method from the 
scale-dependence induced by renormalization group. This amounts to scaling up F(0, yU.) 
from the lower scale /i to, say, mj,. In doing so, we note that the complete solution of the 
renormalization group equation for Oy, 



16/(3/3o) 



+ J2e,i^^)0,i^^) (15) 



with coefficients ej(/i) that describe the admixture of the other operators at 0(a<j), simplifies 
drastically for the matrix elements, since 

{K*^\Oi\B) = forz^r (16) 

at the specific scale and at one-loop level. Thus we have 

«.(m,)^^^/(^^°) 



F(0,mb) =F(0,/i) 



as(/i) 



(17) 



Note, however, that although the initial condition for the solution of the renormalization 
group equation contains seven zeros (the matrix elements in (|TBp), all these matrix elements 
acquire non-zero values in scaling up to and contribute to the rate. In Fig. |l| we show 
F(0,mb), calculated according to Eq. (|r^) and Eq. (|r^) as function of the Borel parameter 
M^. The shaded area shows the working region of the sum rule, 7GeV^ < < 10 GeV^. 
Unfortunately, the sum rule is rather unstable against variation of the Borel parameter which 
is due to a relative sign between the contributions of the quark and the mixed condensate. 
This causes these contributions to be the dominant ones, whereas the perturbative contri- 
bution is rather small (~ 20%), the contribution of the gluon condensate and the four quark 
condensates are negligible (~ 2%). From the figure, we get 

F(0,mb) = 0.75 ±0.10. (18) 

This value is smaller than the one obtained in previous analyses within the QCD sum rules 
approach and larger than the value obtained in [^Q. 

As mentioned before, it would be inconsistent to use this value for the determination 
of the decay rate. The correct procedure is to insert -F(0, fi) as obtained from (|12|) into 
Eq. (P, using the same scale in the coefficient C-j{^). For the combination iVcfepr^ that 



enters the formula for the branching ratio, we use the value 1.81 ■ 10 s obtained in 



which amounts to \Vcb\ = 0.035 for the most recent determination tb = 1.5 ps of the B 



meson lifetime [23|. For the mass of the t-quark that dominates the penguin-amplitude we 
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use the range of values (150 — 200) GeV which is consistent with a recent determination of 
a lower bound on rritop from low-energy data ^4|. Taking all together, Eq. (|) yields the 



branching ratio B{B ^ K*'~f) plotted in Fig. ^ as a function of the Borel parameter M^. 
Again we observe a rather strong dependence on which exceeds all other uncertainties 
in the input parameters of the sum rule, i.e. the quark masses and the values of the quark 
and the mixed condensate. Thus we get a branching ratio with large errors: 

B{B K*-f) = {6.8 ±2 A) -IQ-^ (19) 

which coincides within the errors with the experimental value B{B ^ K*'j) = (4.5 ± 1.5 ± 
0.9). Taking the ratio of the exclusive to the inclusive decay rate, we obtain from Eqs. (^ 
and d): 

^^L_^ = (20±6)%. (20) 

This value again is lower than the QCD sum rule results where values between 28% 
and 39% are given, but it is much larger than old constituent quark model results (4.5% p. 



6% P). It agrees with a recent analysis in p5| where (13 ± 3)% is obtained 



Concluding, we remark that although our results are by far not accurate enough as to 
impose constraints on mtop or possible new physics, we find satisfying agreement of the QCD 
sum rule calculation for B K*'~f with experiment. The main uncertainty of the sum rule 
calculation is caused by different signs in two major contributions which spoils the stability 
of the sum rule in the Borel parameter. Although this uncertainty might be reduced by 
including radiative corrections, their calculation proves prohibitively complicated. 

Acknowledgement: The author thanks M. Misiak for useful discussions about the short- 
distance expansion for b 57. 
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FIGURES 




M_B^2 [GeV^2] 



FIG. 1. The form factor F{0,mi,) as function of the Borel parameter for mi, = 4.8 GeV, 
ms(lGeV) = 0.2 GeV and continuum thresholds s% = 34GeV^, = 1.7 GeV^, evaluated at 
a fixed ratio of the Borel parameters, M^/Mj^^ = 3. The shaded area indicates the "sum rule 
window" . 




FIG. 2. The branching ratio B{B — > K*^) as function of the Borel parameter with the 
same parameters as in Fig. |l]. The shaded area again indicates the "sum rule window" . The solid 
line corresponds to rritop = 150 GeV, the dashed line to mtop = 200 GeV. 
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